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ABSTRACT 


Supersonic flow past two oscillating airfoils is 
analyzed in the thesis using an analytical elementary 
theory valid for low frequencies of oscillation. The 
airfoils may have arbitrary stagger angle. This 
approach generalizes Sauer's solution for a single 
airfoil oscillating at small frequencies in an 
unbounded supersonic flow. 

It is shown that this generalization can provide 
an elementary theory for supersonic flow past two 
oscillating airfoils. This aerodynamic tool will 
facilitate the calculation of pressure distribution and 
consequently the calculation of moment coefficient. 
Torsional flutter boundaries are computed. The results 
for the pitch damping coefficient are the same when 
compared with previous analysis. For arbitrary 
frequencies a linearized method of characteristics was 
outlined. 

The elementary theory that has been developed in 
the thesis can be used for flutter evaluation of 
aircraft carrying external stores. The result of the 
thesis is the derivation of the pitch damping 
coefficient which is necessary to predict the flutter 


Conditions. 
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I. INTRODUCTION 


Aeroelasticity is the study of the effect of aerodynamic forces on elastic bodies. One 
of the interesting problems in aeroelasticity is the stability of a structure in wind. Since for a 
given configuration of the elastic body, the aerodynamic force increases rapidly with the wind 
speed while the elastic stiffness is independent of the wind, there may exist a critical wind 
speed at which the structure becomes unstable. Such instability could create excessive 
deformation and failure of the structure. A particular problem is the FLUTTER of structures, 
which is a self-excited vibration phenomenon. 

Flutter is caused by the interaction of all three forces: aerodynamic, elastic and inertial. 
Consider as an example a cantilever wing mounted in a wind tunnel (with the root ngidly built 
in). Suppose the wing is deliberately deflected and then released. At sufficiently low wind 
speed, the oscillation that follows this disturbance will quickly die out. However, at some 
higher wind speed one observes a steady oscillation, which maintains itself. This is the critical 
flutter speed for the wing. At higher speeds, the oscillation will be rapidly divergent causing 
structural failure. 

When large external bodies such as engine nacelles, fuel tanks, or electronic warfare 
pods are added to the wing of an aircraft, the dynamic charactenstics of an aircraft will be 
changed. Over the last several decades, there has been considerable interest in the calculation 
of the unsteady aerodynamic forces on aircraft. Many of the methods presently in use for 
numerical calculations involve application of a two-dimensional theory to solve three- 
dimensional flow problems. Most of these methods do not have the ability to model an 
arbitrarily shaped aircraft. However, with certain assumptions, it is possible to simplify the 
flow problems to adequately account for the effects of unsteady aerodynamic forces on aircraft 
components. Unfortunately, the transonic flight regime still makes it very difficult to determine 
the unsteady aerodynamic forces because of the nonlinear shock and flow separation effects, 
which occur in this flight regime. 

A good example for the requirement of the wing/stores flutter analysis is the fighter 
aircraft. 


Fighter aircraft are commonly designed for a pnmary mission, such as air supenonty, 
which may require few, if any, wing-mounted external stores. However, for increased 
effectiveness and versatility, many secondary missions evolve which are necessary and require 
the use of a variety of external stores. Thus, many combinations of these external stores must 
be carried at various stations on the wings to achieve the complex, multi-role missions required 
by the operational commands. Sometimes, the wing stores are not available to the aircraft 
manufacturer even during the production of an aircraft and in that case, the flutter analysis 
should be conducted using data not available from the original design. 

A typical modem tactical fighter can carry a great variety of wing stores and 
consequently the total number of possible aircraft/store configurations is huge. Therefore, on 
military aircraft the effect of unsteady aerodynamics is computed for both the clean wings (no 
stores) and the wing with tip missiles. Using these two aerodynamic configurations, the flutter 
analysis is probably done on between 300 and 400 selected wing/store configurations. The 
cost for the above effort would increase by several orders of magnitude if wing/store 
aerodynamics were also considered for each of these configurations. In many cases, this is not 
necessary since the store aerodynamics has a small effect on the flutter speed. However, there 
are cases where neglecting the store aerodynamics will lead to the overestimation of the flutter 
speed. Therefore, the determination of the cost effectiveness of such analysis is a critical step. 

In this thesis, we first review the most important papers, which have been published in 
recent years to analyze wing-store flutter problems. It is evident from the literature survey that 
most approaches are based on modem CFD (Computational Fluid Dynamics) approaches 
which require extensive computations. 

However, there is still some ment in simplified analyses, which reveals the major 
physical effects. Limiting our approach to low supersonic flight speed we propose to replace 
the multiple stores undemeath the wing by a second wing and to apply to a 2-D problem the 
so called “ strip theory approach” which reduces the three dimensional flow problem. This 
simplification has the advantage of accounting for the interference effects between two lifting 
surfaces. Moreover, these interference effects will be stronger than those caused by the actual 
stores. The insights which can be attained with such an analysis, therefore will represent an 


extreme limiting case, but it will nevertheless reveal important physical effects. 


The importance of such interference effects is well known in transonic and supersonic 
wind tunnel testing and in the operation of transonic and supersonic compressors. Platzer 
(1973) presented an analytical theory for the analysis of oscillatory supersonic wind tunnel and 
blade interference effects. This theory is based on the assumpton that the actual airfoil or 
compressor blade can be replaced by a flat plate and that it oscillates with low frequency. This 
makes it possible to expand the solution in powers of the reduced frequency and to retain only 
the zeroth and first order term. Platzer (1973) showed that this theory provides a convenient 
analytical way to estimate the pitch damping as a function of supersonic Mach number and 
pitch-axis (elastic axis) position. 

In this thesis, we modify Platzer’s solution to the case of two airfoils, which are in 
close proximity to each other. Two cases need to be distinguished, namely the one with 
supersonic leading edge locus and the case with subsonic leading edge locus. These two cases 
and the basic theory will be explained in chapter III. The extension of Platzer's theory will be 
given in chapter IV. This is followed by a chapter on the method of characteristic approach, 
which can be used to compare with the analytical results. The final chapter is devoted to the 


flutter analysis. 





II. LITERATURE SURVEY 


EVALUATION OF METHODS FOR PREDICTION AND PREVENTION OF 
WING/STORE FLUTTER. 


Author(s): Pollock, S. J. Sotomayer, W. A.; Huttsell, L. J.; Cooley, D. E 


Source: Collect Tech Pap AIAA ASME ASCE AHS Struct. Dyn. Mater. Conf. 22nd AIAA 
Dya. Spec. Conf., April 6-10 1981, Atlanta, GA, New York, NY, pp. 362-372 


In response to the need to reduce costs and improve safety for flutter evaluation of 
aircraft carrying external stores, the Flight Dynamics Laboratory (FDL) has sponsored several 
efforts in the technical areas of unsteady aerodynamics, flutter prediction, and active flutter 
suppression. This paper discusses each of these three areas as they relate to wing/store flutter 
and presents specific examples from analyses and tests. Steady and unsteady pressure 
measurements were obtained in a wind tunnel at subsonic, transonic and supersonic speeds on 
a fighter wing, tip-mounted launcher and store, and underwing pylon and store. Store flutter 
calculations were performed using both calculated and measured data to determine the 
influence of store aerodynamics on the flutter charactenstics. To improve the accuracy and 
reduce the time and costs of flutter evaluations on the many store configurations carned by Air 
Force fighters, the Flight Dynamics Laboratory (FDL) has sponsored several programs in the 
technical areas of unsteady aerodynamics and flutter predictions. Also, the FDL has been 
sponsonng several efforts to explore the potential of active flutter suppression systems using 
feedback control techniques to provide the required stability and to avoid speed placards. 

This paper reports on some of the FDL research related to wing/store flutter 
prediction and prevention. The research includes an unsteady aerodynamic measurement 
program for a representative fighter wing, with and without tip missile and underwing store 
with test data covering the Mach number range 0.6 to 1.35. Flutter analyses based on an FDL 
computer program specifically for use on aircraft with extemal stores are also described. 
Flutter trends using this computer program are presented for the wing with and without stores 
based on sectional force coefficients from wind tunnel measurements and from theoretical 


calculations. A brief description is given of some FDL programs in active flutter suppression, 
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and typical results are presented for wings with stores, which indicate significant potential for 
improvement in flutter speeds. 

Flutter trends for three different fighter wing/ store configurations were calculated 
using the FACES flutter analysis procedure. Use was made of available modal vibration data 
and measured aerodynamic data The predicted Mach number trends gave minimum flutter 
speeds at transonic and low supersonic speeds as would be expected from the trend of 
measured center of pressure and lift-curve slope data Although the aerodynamics on the tip 
launcher had a somewhat detnmental effect on flutter, the aerodynamics on the tip store had a 
much larger detnmental effect. For the underwing store, the aerodynamic effects were 
beneficial for flutter since the pylon had the effect of decreasing the aerodynamic loading on 


the outer portion of the wing. 


WIND TUNNEL TESTS ON A FIGHTER AIRCRAFT WING/STORE FLUTTER 
SUPPRESSION SYSTEM 


Author(s): Tumer, M. R. 


Source: AGARD, Rep. on a Coop. Program on Active Flutter Suppression, Paper Presented at 
the Structural and Material Panel Meet., 50th April 1980, Athens, Greece AGARD. 


A special system designed using analytical data was tested on a YF- 17 aircraft model in 
the NASA Langley 16' wind tunnel and succeeded in meeting the requirement to increase the 
flutter dynamic pressure by 70% at Mach equals 0. 8. The system was designed using a novel 
procedure, which provides these stability margins, uses minimum control surface movement in 
turbulence, and can be designed using either analytical or empirical data. Two wing tip 
accelerometers and a leading edge control surface were used. Empirical open loop transfer 
functions obtained during the test showed that the analytical data overestumated the response 


of the flutter mode to leading edge control surface excitation. 


AERODYNAMICS USING TIME-LINEARIZATION TRANSONIC FLUTTER 
ANALYSIS 


Author(s): Wong, Y.S., Lee, B-H.K-; Murty, H.S. 
Source: Journal of Aircraft, v. 30, January-February 1993, pp. 144-145 


A survey of the progress made in the development of numencal simulation techniques 
for unsteady transonic flow calculations are presented. Computational methods in three- 
dimensional unsteady transonic flows concentrate mainly on the transonic small disturbance 
equation and time-lineraization approach. An algonthm is introduced for solving flutter 


occurrence. 


WING FLUTTER BOUNDARY PREDICTION USING UNSTEADY EULER 
AERODYNAMIC METHOD 


Author(s): Lee-Rausch, Elizabeth M.; Batina, John T. 


Source: Journal of Aircraft, v. 32, 2 March-Apnl 1995, AIAA, Washington, DC, USA, pp. 416- 
422 


Modifications to an existing three-dimensional implicit, Euler/Reynolds-averaged 
Navier-Stokes code for the aeroelastic analysis of wings are descnbed. These modifications 
include the incorporation of a deforming mesh algorithm and the addition of the structural 
equations of motion for their simultaneous time-integration with the governing flow equations. 
This article gives à brief description of these modifications and presents unsteady calculations 
that check the modifications to the code. Euler flutter results for an isolated 45-deg swept- 
back wing are compared with experimental data for seven freestream Mach numbers that 
define the flutter boundary over a range of Mach number from 0.499 to 1.14. These 
compansons show good agreement in flutter charactenstics for freestream Mach numbers 
below unity. For freestream Mach numbers above unity, the computed aeroelastic results 
predict a premature rise 1n the flutter boundary as compared with the experimental boundary. 
Steady and unsteady contours of surface Mach number and pressure are included to illustrate 
the basic flow characteristics of the tme-marching flutter calculations and to aid in identifying 


possible causes for the increase in the computational flutter boundary. 


APPLICATION OF TRANSONIC SMALL DISTURBANCE THEORY TO THE 
ACTIVE FLEXIBLE WING MODEL 


Author(s): Silva, Walter A.; Bennett, Robert M 


Source: Joumal of Aircraft, v. 32, January-February 1995, AIAA, Washington, DC, USA, pp. 
16-22 


A code, developed at the NASA Langley Research Center, is applied to the active 
flexible wing wind-tunnel model for prediction of transonic aeroelastic behavior. A semispan 
computational model is used for evaluation of symmetric motions, and a full-span model is 
used for evaluation of antisymmetnc motions. Static aeroelastic solutions using the 
computational aeroelasticity program-transonic small disturbance are computed. Dynamic 
(flutter) analyses are then performed as perturbations about the static aeroelastic deformations 
and presented as flutter boundanes in terms of Mach number and dynamic pressure. Flutter 
boundaries that take into account modal refinements, vorticity and entropy corrections, 
antisymmetnc motions, and sensitivity to the modeling of the wingtip ballast stores are also 


presented and compared with expenmental flutter results. 
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WING-STORE FLUTTER ANALYSIS OF AN AIRFOIL IN INCOMPRESSIBLE 
FLOW 


Author(s): Yang, Zhi-Chun Zhao, Ling-Cheng 
Source: Journal of Aircraft, v. 26, 6 June 1989, pp. 583-587 


The flutter of two-dimensional airfoil with external store is analyzed to investigate the 
effects of pylon stiffness on flutter speed. Among the 40 configurations studied, five were 
tested in the wind tunnel to verify the analytical results. The variations of wing-store flutter 
speed si the pylon stiffness can be divided into three types. The curves of the normal and 
flutter frequencies vs pylon stiffness have the same pattern. They can be sketched 
approximately by the aid of the normal frequenaes of the two degenerated two-degree-of- 
freedom systems, i.e., and the freely hinged and nigidly connected store cases. A limiting flutter 
speed for very small pylon stiffness is deduced, which is useful to identify which type of flutter 
the configuration studied belongs to. 
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STUDY OF THE EFFECT OF STORE AERODYNAMICS ON WING/STORE 
FLUTTER. 


Author(s): Turner, C. D. 


Source: Collect Tech. Pap. ALAA ASME ASCE AHS Struct Struct. Dyn. Mater. Conf. 22nd 
ALAA Dyn. Spec. Conf., Apal 6-10, 1981, and 1981, Atlanta, GA, AIAA (CP 812), New York, 
NY, pp. 343-351, Paper: 81-0604 


This study represents the first systematic analytical study of the effect of store 
aerodynamics on wing/store flutter. A large number of wing/store single carnage 
configurations and parameters were ease in the study; multivanate analysis techniques 
were used for the first time to analyze wing/store configurations, modal data, and flutter 
results. The results of the multivariate analysis indicate that it may not be possible to develop 
general guidelines, but it is possible to develop specific guidelines for use with a particular 
aircraft. This study was the first attempt to do a systematic analytical study of the effect of 
store aerodynamics on wing/store flutter. To determine this effect flutter analyses were done 
on four aircraft with single carnage of three basic store types. In all 308 configurations were 
analyzed with and without store aerodynamics (tip missiles, tip tanks, and underwing stores). 
The effects of stores have been analyzed. The results of the factor analysis indicate that it may 
not be possible to develop general guidelines, but it is possible to develop specific guidelines 
for use with a particular aircraft. The conclusions of the study, as far as it concems the effect 
of store aerodynamics on wing-store flutter, show that 60% of the tested configurations 
require a change in the flutter speed no more than 7%. Therefore, there is no need to evaluate 
the requirements for aerodynamic modeling of the store. However, for the rest of the 
configurations there is such a need. Moreover, 75% of the last reveal that nonrealistic flutter 
results are obtained without taking into consideration the store aerodynamic. Of course, the 


results of this study represent the configurations that were used to generate the data base. 
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AERODYNAMIC MODELING OF AN OSCILLATING WING WITH EXTERNAL 
STORES 


Author(s):, Sotomayer W. A., Dusto A. R, Epton M.A., Johnson F. T. 


Source: AIAA , New York, NY, pp. 245-252, Paper: 81-0609 


An analysis of the steady and unsteady aerodynamic forces acting on a fighter aircraft 
wing with stores has been done. Computations were performed with paneling methods 
capable of presenting arbitrary aircraft configurations in subsonic and supersonic flow. 
Interference effects from a tip store and an underwing pylon store in varying stages of 
completion were also analyzed. Detailed comparisons between experimental data and 
numerical computations are also made. 

In order to calculate the pressures, forces and moments on an aircraft a numerical 
approach is proposed. Numerical computations are based on solving an integral equation 
formulation of the flow problem being considered. Steady state calculations were done with 
two numerical methods. The first of these 1s a pilot code developed by Johns. In this method, 
distributions of linearly varying sources and quadratically varying doublets are used to 
represent the aerodynamic surfaces. Wakes are represented by doublets with constant 
streamwise strength and linear variation in the spanwise direction. A method developed by 
Woodward was also employed in making numerical calculations. This method utilizes 
distributions of sources and vortices to represent an aerodynamic surface and the wakes shed 
from various components. 

For unsteady flow, each of the methods makes use of sources and doublets as the basic 
aerodynamic singularities. In a method developed by Johnson for unsteady subsonic flow, 
distributions of sources and doublets are used to represent the aerodynamic surfaces. In the 
wake, doublets of fixed strength and position are used to represent the effects of unsteady 
wake motion. Calculations were also done with a special doublet lattice method. In this 


method, a wing is represented as a sheet of doublets, components such as the fuselage or a 
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store or nacelle can be represented with an axisymmetric distribution of sources with doublets 
on the surface. 

Approximations are imposed on these distributions and on the kernel functions used to 
represent their influences on the flow. These approximations are briefly described along with 
the flow boundary conditions imposed on the aerodynamic surfaces and wakes. 

It is interesting that the small disturbance partial differential equation is being used for 
unsteady inviscid incompressible flow along with the linearized form expressed in terms of the 


velocity potential @. In addition, simple harmonic motion is assumed as in the thesis. 
p=Refp*e®] where @* is a complex quantty. 


Moreover using the Helmoltz equation a more simple relationship is established for @* with 
the form of an integral equation. An approximate solution of the last equation permits an 
evaluation of the entire flowfeld. Numerical solutions of the last equation are achieved 
through an aerodynamic influence coefficient method. The boundary conditions are 
established taking into consideration the configuration of the F-5 aircraft, a relatively old 
aircraft. The aircraft is comprised of the following components: 

1) wing 

2) underwing pylon 

3) missile rack 

4) missile body. 

As in the thesis, the above surfaces are approximated as infinitesimally thin surfaces so 
the thin wing theory applies except the for missile body. The boundary conditons are 
expressed in terms of Taylor series expansion about reference surfaces taking into 
consideration the thickness. 

Initially boundary conditions for steady state flow are made use of and unsteady 
boundary conditions are taken as a perturbation about a steady state condition. In addition, 
variations in Mach number, frequency of oscillations, and interference effects arising from 


component build up are presented. 
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For the wing missile, the boundary conditions are of mixed type (Dirichlet and 
Newman). On the surface of the missile, the boundary conditions are similar to the other parts 
of the aircraft. 

Using a Woodward method, a lifting surface is divided into numbers of aerodynamic 
panels number, containing distributions of sources and vortices. A source distnbution 
represents a fuselage [or a pod or external store]. Wing thickness 1s represented by a linearly 
varying source distribution in which the strength is equated to the chordwise slope of the wing 
thickness. Camber, twist, and lifting effects are represented by a linearly varying vortex 
distribution where the strength is determined to satisfy tangential flow at panel control points. 
An iterative procedure is employed in solving the boundary value problem. For analysis 
problems, the surface slope is described, and singularity strengths are determined by inverting 
the matnx of aerodynamic influence coefficients. With the strengths of the aerodynamic 
singulanties known, the u, v, and w velocity components at a given point may then be 
determined. Pressures, forces, and moments are calculated by numerical integration. 

In the doublet lattice method, the aerodynamic surface is subdivided into a senes of 
infinitesimally thin panels. Along the quarter chord line of each of these panels is contained a. 
distnbution of acceleration potential doublets. The strength of the distribution is constant but 
is not known. Specification of the normal velocity at a set of points on the surface determines 
the loading on each element. Locating the lifting elements at the quarter chord and the 
collocation point of the three-quarter chord of the midspan of each panel usually results in 
reasonable success. An approximation for the lifting pressure coefficients in terms of the 
induced normal velocities may then be found. Lift, moment, roll and generalized forces may | 
then be calculated. 

Recently, a series of wind tunnel tests were conducted by NLR of the Netherlands with 
sponsorship from the Air Force. These tests involved measurements of steady and unsteady 
pressures and forces on a model of the F-5 wing with external stores. Several external store 
arrangements, which are described and shown in this paper, are compared with numerical 
calculations. Experimental results and numerical calculations show that: 

1) For subsonic flow, the interference due to the tip store is much greater on the 


outer portion of the wing than it is on the inner portions. 
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2) The interference effect due to the pylon store is significant on the lower surface of 


the wing for subsonic flow; there are also noticeable effects on the upper surface of the wing at 


subsonic speeds. 


3) For supersonic speeds, the zone of influence remains localized within the Mach 


cone. 


In general, numerical calculations and expenmental results agreed reasonably well with 


each other. 
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HI. PROBLEM FORMULATION 


The interference between the airfoil (wing airfoil) and the wing stores is examined by 
making certain assumptions that simplify the problem. The wing airfoil is replaced by an 
oscillating flat plate airfoil and the wing stores by another oscillating flat plate in distance d 


from the wing airfoil. Two cases need to be distinguished, namely: 


SUPERSONIC FLOW PAST TWO 
E OSCILLATING AIRFOIIS 


WING AIRFOIL 


II E 





II* 


STORE AIRFOIL x 


Fig. 1 


1) Supersonic flow past two oscillating airfoils with supersonic leading edge locus 
shown in Fig. 1 (above) 

2) Supersonic flow past two oscillating airfoils with subsonic leading edge locus 
shown in Fig, 2 
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Fig. 2 


The differences are obvious and will be further analyzed dunng the examination of 
each case. In both cases multiple reflections occur. As explained 1n the preface the interference 
between the two airfoils under certain conditions could result in unsteady flow phenomena like 
flutter. 

Therefore, the study of the oscillating airfoils is important in predicting the flutter 
characteristics or dynamic response of the wing airfoil when different stores are attached. For 
the case shown in Fig. 1: 
tan f « cote (III-1a) 

Where the case shown in Fig. 2 requires: 
tan J > cota (III-1b) 
The flow 1s assumed a non-viscous compressible two-dimensional flow of a perfect gas 


govemed by the continuity equation, the Euler, and the energy equation. The continuity 


equation 1s: 

A pee © | 5 (1-2) 
Dt Ox Oy 

The Euler Equations: 

eu + oP, =) (HI-3) 
Dt po 

Dv 1, (111-4) 
Dt poy 

The Energy Equation: 

Ds 

—=0 II-5 
Di (11-5) 
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Because of the assumption of small amplitude oscillations of the airfoils, all flow quantities are 
considered small perturbations linearly superimposed on free stream quantities. The velocities 


can therefore be wntten as: 


u =u. +u (TII-6) 
C=C, +c (11-7) 
ver (TII-8) 


In addition, the pressure and density perturbations are linearly supenmposed on freestream 


quantities. Therefore: 


s pp. (I1I-9) 


ADDE (111-10) 


The local velocity of sound is given by: 


= 2 
£ 2) RES (III-11) 
Furthermore: 
£ = const (11-12) 
p 
Taking the total differential: 
l ] 
Or: 
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= yf (III-14) 
dp P 
So: 
Pog” (III-15) 
p 


From AA 4318 class notes, the surface boundary conditions will be calculated. Requirement 1s 
that no flow will penetrate the solid surface of the airfoil resulting in the flow tangency 


condition: 
V-7.)-n=0 | (1-16) 


where V,V, are the velocities of the free stream and the surface over which the fluid flows.If 


the body is described by the equation: 
F(x, y,t)=0 (11-17) 
then the normal to the surface at any point is given by the following relation: 


VF 


"M 


(11-18) 


Taking into consideration that the substantial derivative of a surface particle should be zero: 


AV, VE =0 11-19) 


one obtains from (11-16): 
VNF=V,-VF (III-20) 
and therefore: 


2 


DE ae (I-21) 
D ôt 


Or for two dimensional flow: 


ici En o + yo = 0 (III-22) 
Dt à Ox Oy 


Particles in contact with the surface must have the same normal velocities as the surface. This 


means that the rate of change of F is zero hence: 


DF 


= 11-23 
p (111-25) 


For an airfoil the equation of the upper surface can be written as follow: 


F(x, yt) = y— yy (x,t) = 0 (111-24) 


Where y, is the distance from the chord line to the upper surface. The equation of the lower 


surface can be written 1n a similar way: 


F, (x,y,t) = y — y, (x,t) = 0 (111-25) 
At y=y,: 


Dri La E ma — (III-26) 
Dt 
A 


DH, __ = cL s I 
Di -(u) +v=0 (III-27) 


Yu Oy 
RU II-28 
ay oy (111-28) 
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Thus, the normal flow velocity can be written: 


nn (I-29) 
Ot Ox 


Applying the assumption of linear perturbation theory one has u = u +u , therefore: 


Oy, u Du u 


yv = — —=. AX 11-30 

_ y = y, (x,t) (11-30) 

PER CI, = y = (xt) (111-31) 
ot Ox 


Neglecting the higher order terms one obtains: 





Oy Oy 
= => E y= if II-32 
v= = tu. 2, y= y, (x, t) (III-33) 


This normal velocity will be expanded as Taylor series of the normal flow velocity around 
y=0: 


v(x t)= v(x,0* ,t)+ y, >t y Ploe), e II-34 
ME 21 Ə2y (I ) 


v(x t)= v(x,0* E y PELO, >L y (0 1), iun 11-35 
JL» 21 O” y (I ) 


Using the assumption of thin airfoil and small linear perturbation theory the higher order 


terms may be neglected and one obtains: 


v(x, y,,t)= v(x,0*,r) (III-36) 


v(x, vn ,t) = v(x,07 y t) š 


(11-37) 

Therefore: 

0 = Dern Du y= 0" (11-38) 
ot Ox 

x,0°)= Peay, Me y=0 (11-39) 
ot Ox 

For simple harmonic motion of the airfoil: 

y, 7 he (11-40) 

ea hx e (HI-41) 


Therefore, the lineanzed flow tangency equations can be written: 


v(x,0*,:) = [i@h,(x)+u, - Jg ^ URN (I1I-42) 





v(x,0",2) = [ioy, (x) +u, Pee, vl ` (HII-43) 
x 


The above relationships describe the normal velocity on the lower airfoil (store airfoil). 


Similarly, one obtains for the upper airfoil (wing airfoil) at y=d: 
vx, d - it) = [i@h, (x) +u, Aye y=d (11-44) 
x 


if the store airfoil executes a pitch oscillation about x = x, one has: 


y 7 -0(x - x, e^ (TII-45) 


v(x,0*,£) 2 [74,0 — io0(x — x, Je" at y - 0* (III-46) 
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and similarly for the wing airfoul: 
E (11-47) 
v(x,d”,1) - [- Olu, cosó — @.sin ó(x — B — x} ¡Ola sinó + 0w.cosó(x-B- Da 


(11-48) 
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IV. ELEMENTARY THEORY 


The case of supersonic leading edge locus will be examined assuming first that both 
airfoils are oscillating slowly in a supersonic flow. The airfoils execute vibrations of identical 
modes and amplitude but with different phase angle 6 between them. It is obvious that the 
flow field is such that no disturbance can propagate upstream of the airfoil leading edges. 


Consequently, the perturbation potential can be wntten as: 


D(x y,t) =0(x,y,k) Eu (IV-1) 


where x, y are non dimensional coordinates and t is non dimensional time. The differential 


equation for the perturbation potential in a non-dimensional form then 1s: 
coto-0, —@, +2ikM°p,—k *M°*p=0 (IV-2) 

The hnearized boundary conditions are: 

9,&0)=v@&) on the lower(store airfoıl) at y=0 (TV-3) 
where v(x) 1s the downwash amplitude. 

DRS ze Sy(x-dtanf) on the upper(wing) airfoil at y=d (IV-4) 

where it is assumed that the two airfoils oscillate with a phase angle 6. 

Assuming oscillation around x, v(x) becomes: 

v(x)=O[1 +ik(x—x,)] and without loss of generality we can set 071 (IV-5) 


For small oscillations (low frequency) the potential can be expanded and neglecting the 


higher order terms of k one can write: 


Q(xy,) 2X(xy) Fk (xy) | (IV-6) 
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Therefore the differential equation (II-2) becomes: 
cot a-X_— X,=0 (IV-7) 
cota P, —P_+ 21M"X,=0 (IV-8) 


A general solution of the above system of equations has been given by Sauer (1950). It reads: 


X(xy)=8(2) (IV-9) 
‘YP (x,y)=h(z)+iMyeg(z)/cosa where z=x—ycota (IV-10) 
T(xy)- h(z) +iMy g(z)/ cosa where z=x+ ycota (IV-11) 


h(z) and g(z) are arbitrary functions for positive arguments of z and zero for z<0, it is obvious 


that the expressions of the above functions should be consistent with the last assumption. 


Using the two solutions for left and nght running Mach waves in the supersonic flow 


field of slowly oscillating airfoils pressure coefficients can be obtained. In order to do so the 


flow field between the two airfoils has to be divided into several zones. The number of zones 


depends upon A. It is obvious that for A>1 where A=dcota there is one zone along each 


airfoil, for 0.5<A<1 there are two, for 0.33<A<0.5 there are three and so on. In the first zone 


I there is no interference from the airfoil wing (upper wing), which means that the lower airfoil 


(zone I only) does not sense the upper airfoil. In that case, the solution has the form: 
Puy k)=8(2) + {h(z)- iMyg(z)/cosQ} 

The boundary conditions that must be satisfied at y=0 are: 

A,=—cotal-g (z)=-1 

P,=—cotah'(z)1Myg(z)/cosa==(x —xy ) 


From the last two equations it is obvious that: 
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g(z)=ztana 

h(z)=—iztana(x)+ztan7a/2) 

Since we are assuming a slowly oscillating airfoil, the solution to the PDE has the form: 
«oy k)-g(z) kh (z)-iMyg(z)/coso] 

where z=x-ycota 


In order to calculate the pressure distnbution on the upper surface of the lower airfoil in the 
zone I the following relation should be used: 


cp=-2[ıkp +9, ] 

Substituting the known functions @ and ç, 

p=ztana+k {-iztana(x,+ztan”/2) —(Myztana/cosa)} 
0,=d/dx[ztana +k{-iztana(x,+ztan 0/2) -GMyztana/cosa)}]= 
=tana—iktana[x,+ztan“a.—(Myz/cosa)]|= 
=tana—ik[tana(x,+ztan“a)+( y/sinacosa) (sinaL/ cosa)]= 
=tana—ik[tana(x,+ ztan ct) + (y/ cos 0)] x 

=tana-ikftana (x,+ (x-ycoto)tan a+ (y/cos*a)]= 

—tanat-ik [tana (x, *x tan A-ytana+ (y/ cos'a)]- 
=tana-ikftana(x,+ x tan a-ycota]- 


Therefore, the pressure distribution is given by the following relationship: 
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cp7-2(tana-ik [tana (x, *- xtan^o—yco ta} t+ik[ztanatk {—iztana(x,+ ztan 0/2)-¡Myztano/ 


cosa}]} 


Taking into consideration that higher order terms of k are neglected the pressure distribution is 


given by the following relationship: 

| cp=2tana{—1+ik|x,+x (tan’a—1)+2ycota]} 
and therefore on the airfoil for y=0: 
cp=2tana{-1 +ikfxotx (tan’a-1)]} 


In the zone I* the perturbation potential can also be waitten directly from Sauer’s single airfoil 


solution: 

Py.) =gr(21)+k[h,(21)+ ¡Myg(21)/cosa] 

Where: 

z¡=x— Á—B +ycota 

The boundary conditions that must be satisfied at y=d are: 

X =cota.g’(z) = E 

V, zcotochy'()-Mgi (n)/cosoi My gi' (6) 736 (s--B): 
After integration and algebraic manipulations: 

gi (z,)=~z,tana-e® 

hy (2,)=iz,tanc-e” (x +2, tan2o/2+Md /cosa) 
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where again the condition has been imposed that the functions g;(z;) and h,(z;) are zero for 
negative and zero arguments. 
Following the same procedure as in zone I, the perturbation potential @ will be 


calculated substituting the known functions g;(z,) and h,(z;): 


P=—2; tane ik (tana Eo nc 2+Md/ Mae) cosa} 


9,=-tane-Hiktancı e { (xp+z,tan°at+Md/ cosa)-My/cosa}= 
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=-tana-e+iktana-e® [(x+(x-A—B+ycota)tana+Md/cosa)-My/cosa)= 


scan oce ( (xo +(=—B)tana—A tan a+ycota+d/sina-cosa) 


=—-tana e 
-y/sinacosa} = 


eee OL: . { (xo+(x—B) tan oct A-ycota} 


=-tana-e 


Because d=A/cota , tana—1/cosa-sina=cota 


This gives the following result for the pressure distnbutions on the lower surface of the upper 
airfoil (wing airfoil) in zone I*: 


Ô 


cp=2tana-e’ —2iktana-e { (x)+ (x—B)(tan’a-1)+2A-—2ycota} = 


cp=2tanare{1-ik { (xp (x-B)(tan’or-1)+2A—2ycoter} 


and therefore at y=d: 


cp=2tanarre! {1a { (xo (xB) (tan2a-1)} 
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Assuming that 0.5<A<1 in addition to zones I, I* there are two more zones II, II*. In 
zone II the perturbation potential is due to the initial waves from zones I, I* plus the reflected 
wave at zone ILSo the reflected waves must be determined by adding the functions 9(z>) and 
h»(z;). In the zone II Sauers solution for a single airfoil applies and following the same 


procedure as previously: 

9.5) =g(2)+k[h(z)-iMG8(2/cosot ] +g (z0+k[hi(z0+iM(g (z)/coso]+ 
82(Z2)tk[ho(Z2)-iMygo(z2 )/cosa] 

The variable z has the following form: 

z2=x-2A+ycoto 

This ensures that z,<0 upstream of the leading edge of the airfoil. 


The boundary conditions that must be satisfied at y=d are: 
x — , f f EN ¡0 

y" cotacg' (z)-cota-g;' (zj)-cotac-g»' (z2)7—e 
'P,—7-cota-h;'(z)*cota-h;' (z;)-cotac-h? (z2-M [g(z)-g1(z1)-g(z2)]/coso-iMdcota [g: (z1) 
+g1 (1) +02 (22)]/ coso= ie (xx) B) 


After integration and algebraic manipulations similar to the previous calculations for zones I, 
FE 


gol Z2)=ztana 
ha(z2) =— izotana Got Z2/ 2tan^a +2Md/ cos) 


Therefore, the perturbation potential is given by the following relationship: 
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Ô 


,=2tana—tana-e’ —iktana(xp+ zotan a+2Md/co SQ) "E ‚tan’a+Md/ cosa) 


~iktana(x,+ ztan*o )-ıkMy[tana-tana+ tana-e® ]/cosa= 
—2tana-tano-e P iktana[(9* z;tan o: 2Md/ De ztan a +Md/cosa)-(x,+ ztan OL 
yeMye /cosa]= 


Ô 


=2tana—tana-e—iktan a[x,+(x-2A+ycota) tan’a+2Md/ en (x-A-B+ycota) tana 


* Md/coso)-(x,* ( x-ycota) tan o) +Mye?/ cosa]= 
=2tana=tagore O + (x-2A) tan o. 2Atana/ A e (x—A—B+yco tot) tan “cr 
+Atana/sina-cosa)-(x +x m sinQ-cosa)]= 


—-2tano.-tano-e  iktano[x, -xtan a 2A -e (xt (x-B+ycota) ES 


/sina-cosa]= 


=2tana-tana-e°-iktana[2x,+2xtandat2A-edx,+ (x-B-ycota)tan^o-- A) ye" /sina-coso]- 


Ô 


—2tana-tana-e —iktano[(2-e9 Y(x,xtan oc- A) s (Bus as ycota)|= 


-tano.(2-e ^ ik [Q9 onse EN Mac ycoto)]j 


This gives the following result for the pressure distributions on the lower surface of the 
upper airfoil (wing airfoil) in zone II: 


Cp--tana( DR (x;+(&-B)(tan’a-1)+2A-2ycota)- 


2(xyx(tan a—1)+2A)) 


and for y=d we have: 


Cp=-tana(2-e +ikfe + (x-B)(tan^a.-1))-2 (x,* x (tan a.-1)-2A)) 


In zone II* the perturbation potential is due to the initial waves from zones I, I* plus 
the reflected wave at zone II, II*.So the reflected waves must be determined by adding the 


functions g3(z;) and h3(z;). Therefore, in the zone II* Sauer’s solution for a single airfoil 


applies and following the same procedure as previously: 

PEYE) =g(2)+k[h(2)¡MGeg(2)/cosa ]|Tgi G) +kfhy (21) +My, (1)/coso] + 

83(25) Fk[hs(3)-3Mygs(z; )/ cosa] 

The variable z; has the following form: 

23=x-A-B-ycota 

The boundary conditions that must be satisfied at y=0 are: 

X, =— cota-g'(z) *coto-g;'( z;) -cotocga' ( zj)* -1 

‘Yy=—cota-h' (z)+cota.-h,’(z;)—cota-hs' (z3) AM[g(2)7g: (21) * g3(23)]/ cost -— i(x—xo) 


After integration and algebraic manipulations similar to the previous calculations for zones I, 
I 


g3( Z3)=—z,tana.-e” 
h3(z3) =iz3tana, (xp+z3tan"a/2+Md/cosa)e” 


Therefore, the perturbation potential is given by the following relationship: 


Ô 


P,=tano-2tanct-e +iktana (fx +2, tano +Ma/ cosa) te” Got z;tan a+Md/co so) (Xp + 


ztan" o. )-My/cosa} 
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= iktana e [xy + (x—A—B+ycota) tan a +Md/ cosa] À [xk 


(),=tana—2tana-e’ 


(«-A-B-ycota)tan’a +Md/cosa] -(x,+ (x-ycota)tan’a )-My/cosa} = 


Ô 


Q,—tana.-2tana.e' ik (eP2tano [xo (x—A-—B)tan^a-- Md/coso]-tana(x,- (x-ycota)tan a)- 


y/cos’a} 


Ô 


p,=tano-2tanc-e"”+ik (rana [xo + @-B) tan’a-Atan’a+ Adtana/sina-cosa]-tana(x,t 


(x-ycota)tan^a )-y/ cos’a} = 


Desi (ePotano px (x—B) tan ot A]-tana(x,* (x-ycoto)tan a)- 


P,=tano-2tanc-e 


y/cos’a} = 


9,=tanct-2tancr-e+ik {e 2tanctfxo+ (xB) tan"o-+ A]—tance(x + xtan at ycota)} 


Because: y/cos’a+ ycotatan a.——y 


In order to calculate the pressure distribution on the upper surface of the lower airfoil 1n the 
zone 11* the following relation should be used: 


cp =-2[tkQ +0, ] 
Substituting the known functions (9 and 9, 


=2tano—Atanot- e —2ik {e®2tanafxo+(x-B)tan’o+ A]-tanaı(&,+xtan?a-+ ycoto)} 


—2ik{ztana—z,tana- Eos ZjtanQat- eo higher order of k terms which will be neglected 
Therefore: 


O D Oy n (x-B)tan^a-- A]-tana(x,xtan a+ ycota)} 
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Cp=2tana—4tana-e! 


i iÓ 
—2ik { (x-ycota)tana—( x—A—B+ycota)tana: e 07 x—A—B-ycotQ)tana ey 


After several relatively simple algebraic manipulations: 


Cp=-2tanaf 1-2e°-ik{2efx,+&-B)(tana-1)+2A]-&,+x(tan?o-1)+2ycota)}] 


Since the pressures have been calculated for the flow field between the airfoils a similar 
approach will be utilized to obtain the pressures on the bottom of the lower airfoil. In that 


case, the solution has the form: 


Q(x y,k) = g(z) +k{h(z) tiMyg(z)/cosa} 
The boundary conditions that must be satisfied at y=0 are: 


X, =cota-g'(z)= —1 

‘Y= cota-h’(z)+iMg(z)/ zn —,) aty=0 
From the last two equations it is obvious that: 
g(z)=ztana 

h(z)=-iztana(x+ztan?a/2) 

The solution again has the form: 

Gy) -g) Fk[h(z)* iMyg()/cosa ] 


where z = xtycota 


In order to calculate the pressure distribution on the upper surface of the lower airfoil in the 


zone I the following relation should be used: 
cp ——2[iko *9,] 
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p=-ztana +kf—iztano(x+ztana/2)+(¡Myztana/cosa)) 
p,=d/dx[-ztana.+k(—iztana(x+ztan"0/2)+(¡Myztana/cosar))]= 
=—tano— iktana[x+ztan a—(Myz/cosa)]= 
=-tana-ik[tana(x,+ztan’a)—{y/ sinacosa)(sinat/cosa)]= 

=tana — 1k[tana(x ¿+ ztan a) +(y/ cos’a)]= 

Therefore, the pressure distribution 1s given by the following relationship: 
cp=2tana{ 1—ik[tana(x)+ x (tan a-1)+2ycota} 


Using the results for zones I,I*,II,II* it 1s possible to calculate the pressure jump across 


the bottom airfoil (wing store airfoil) from the following relationship: 


Acp(x,0)=cp(x,0°)- cp(«,0) 


Two integration intervals, shown in Fig.4, can be distinguished in the lower airfoil taking into 
consideration that 0.5XAX1: 

1) O<x<A+B 

P A b= A 


In the first integration interval: 

Acpı(8,0)=cp(x,0”)- cp(&,0 )= 

=2tanat l-ik[tana(x,+x (tana—1)+2ycota) 

+2tano4 1-ik[tana(x,+x (tan“a—1)+2ycota) and y=0 
Acp,(x,0)=4tana{ 1—ik(x,+ x(tan*a-1)} 


In the second integration interval: 
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Acpy(x,0)= 2tano.{ 1—ik[tana(x,+ x(tan°a-1)+2ycota} 


*2tano[1-269 ik (269 «4 GB) (an^ 0.71) 2A]- (x, x (tan 0.1) 2ycoto)] 


=+2tanof1—26-+ik (26 fxo+ («-B)(tan’a-1)+2A]-2(&;+x(tan’a-1))}] since y=0 


But e”=cosö+isind and the last expression can be written: 


-2tana[1-2(cos8*isinB) ik {2 (cos6+isind)[xo+ (xB) (tan"a—-1)+2A]-2(x+x(tan“o-1))}] 


This expression is a complex equation with real and imaginary part. In order to calculate the 


pitch-damping coefficient only the imaginary part of Acp is needed. Therefore one has: 


ImAcp,(x,0)=-4tana (x, +x(tan’o-1)) 
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ImAcp, (x,0)=+2tana[-2sind/k+2cosö[xo+ (x-B)(tan’a-1)+2A]-2(x,+x(tan’a-1))] 
The lift and moment coefficient can be calculated from the following relationship: 
C; —ol 'Acpdx 

Cal “Acp(x—x,)dx 

By definition: 

C,=9[ C,e ik Cue] 

Ge Curt Cue 

(1/ik) f**PImAcp,(x,0) (x-x,)dxt asp! "ImAcpg(x,0) &—x;)dx- 


of *P_-4tana(xy+x(tanO-1)(x-x,)dx+aspl2tano[—2sin8/k-+2cos8{xo+ (x—B)(tan°a—1)+2A]—2( 


xy x(tan^a.1)) (x-xydx 
The last two integrals will be calculated separately. The first one is: 
2tana[-2x,? (A+B)+ x,(A+B)(A+B)-2(tan’a-1)[x,(A+B)?/2-(A+B)’/3)] — (IV-12) 


The second one can be calculated using Maple. The results are a complicate expression. 


Two calculations have been included for the convenience of the reader since in both integrals 


there is a common factor 2tana. The result is: 


2tana {2x,+tan’a-B’A?B-cosd-tan?a-B’/3+2tan’a/3-sind-A?/k-sind-B?/kt+sind/k+2cosd/3 
+2x,sind-A/k-2x,’cosd-B-2x ?cosd-A+4x,cosd-A+2x,cosö-B+2x,sind-B/k+2x, A+2x, B 
+2cos&-tan’aA?/3+2cosd-AB?-4x,AB-2sind-AB/k-2/3+4cosd-A?/3+2B?/3+2AB’+2BA? 
-2tan’a-AB°-2tan’a- AB?-2cosd-tan?a,/3-2x,cosd-x,tan?a-cosd-B-2cosd-A-2x,sind/k+2x,” 


sind/k+2x, cos6-2x, +x,tan’a-A”+x,cosö-tan’atx;tan?a-B’+cosö-Btan’a-2x,cosö-Btan’a+ 
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3cos6-A7B—2tan’a-A’/3—2tan*a-B’/3+cos6-B’/3+2A’/3—2x,B°+2x,tan*a-AB 
—2x,A’—2x,cos6:BA+x,cos6-tan’a@-B*—x,cos5-A*—x,cos6-tan*a-A*} (IV-13) 


Therefore: 


Ce =2tanaf-2x, (A+B)+ x,(A+B)(A+B)-2(tan’a-1)[x,(A+B)?/2-(A+B)’/3)]+ 

2tana {2x,+tan’a-B’A’B-cosö-tan’a-B°/3+2tan?a/ 3-sind-A?/k-sind-B*/k+sind/k+2cos6/3 
+2x, sind-A /k-2x,cosö-B-2x, cosö-A+4x,cos6-A+2x,cosö-B+2x,sinö-B/k+2x, A+2x, B 
t2cos0-tan a.A?/34-2cos6-AB/—4x,A B —2sin6-AB/k—2/3--4cos6-A? /3- 2B? /5-2AB^-2BA? 
-2tan’a-AB°-2tan’a-AB’-2cosö-tan”a /3-2x,cosö-x,tan”a—-cosd-B-2cosd-A-2x,sinö/k+2x,” 
sindö/k+2x,’cosö-2x, +x,tan’a-A”+x,cosö-tan’a+x,tan”a-B”+cosö-Btan’a-2x,cosö-Btan’a+ 
3cosö-A”B-2tan’a- A? /3-2tan’a-B’/3+cos6-B’/3+2A°/3-2x,B”+2x,tan’a-AB 
—2x,A°—2x,cos6-BA+x,cos6-tan°aB*—x,cos6-A*—x,cos6-tan*a:A’} (IV-142) 


This is the major result of our investigation. It represents the pitch damping coefficient of the 
lower airfoil in the presence of the upper airfoil. Note that both airfoils are oscillating in pitch. 
Also the two airfoil oscillation may lag by a phase angle 6. 

If B=0 and 5=180° the above result for the pitch damping coefficient is the same as 
that of an oscillating airfoil mounted close to a stationary airfoil or wall at a distance d/2(M. F. 
Platzer, H. G. Chalkey). As will be shown later the same result can be derived using the 


elementary theory. The pitch-damping coefficient for that case 1s given by the relation: 
Cye-——2tana (x, (4—2A)—x,(3A^—4A -2-- (A?—2)(tan^a.—1))--2A)-2A —2(2—A?)/3 (tan^a.-1)] 
(IV-14b) 


Where A=2hcota, and h is the airfoil distance. 
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Fig. 5 


For th ! 
u € m case of zero stagger B=0 and out-of-phase oscillation d=180 the 
general result equations (IV-12) and (IV-13) simplifies to (IV-14). 
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This agrees with the result previously obtained by Platzer and Chalkey (1972) if 
A=d.cota in replaced by A=2hcota in replaced by A=2hcota used by Platzer and Chalkey. 
Two unstaggered airfoils oscillating out-of-phase are equivalent to a single mounted at a 
distance h=d/2 from a solid wall of a non-moving airfoil. 

It can be seen from Fig. 5 that the interference is highly destabilizing with increasing 
amounts of interference. 

Note that A=0.8 implies a larger amount of interference than A=0.9, for example. At A=0.9 
instability is encountered from low supersonic flow to a maximum of M=3, which at A=0.8 


instability occurs up to much larger supersonic Mach numbers. 





DAD AA Do 08.2 4 


Fig. 6 
Another interesting case occurs for in-phase oscillation 6-0 of two unstaggered airfoils. The 
instability boundaries are shown in Fig. 6. Again, it is seen that increasing amount of 


interference leads to increasing regions of instability. 
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Fig. 6 shows curves of zero pitch damping coefficient (with torsional stability 
boundaries) for an airfoil mounted close to another airfoil when B=0 but oscillating in phase. 

Following the same procedure a special case of the above will be calculated. In that 
case the upper airfoil is assumed to be stationary while the lower airfoil (store airfoil) oscillates. 


The results will be compared with previous analyses (M.F. Platzer, 1971). 


Starting again with Sauers solution for the slowly osallating airfoil in unbounded 
supersonic flow it is shown how this solution can be extended to consider the upper airfoil 


(wing) interference. To further simplify the comparison it is assumed that: 
JM? -125 « 
E 


Therefore, the moving airfoil is mounted at half distance from the stationary upper airfoil in 


compare to the previous analysis. The solution 1s given by the following set of equations: 
X&y)=g(Z) 

(x,y) =h(z)+iMyg(z)/cosa where z=x—ycota 

‘Y= h(z)+iMy g(z)/cosa where z=x+ycotal 


h(z) and g(z) are arbitrary functions for positive arguments of z and zero for z<0, it is obvious 
that the expressions of the above functions should be consistent with the last assumption. 

Using the two solutions for left and nght running Mach waves in the supersonic flow 
field of slowly oscillating airfoils pressure coefficients can be obtained. In order to do so the 
flow field between the two airfoils should be divided into several zones. The number of zones 
depends upon A. It is obvious that for A>1, there is one zone along the airfoil, for 0.5<A<1 
there are two. So there are three regions. In the first zone I there is no interference from the 
airfoil wing (upper wing), which means that the lower airfoil (zone I only) does not sense the 
upper airfoil In that case, the solution has the form: 


Py) = gz) +k{h(z)- iMyg(z)/cosa} 


The boundary conditions that must be satisfied at y=0 are: 
X, =—cota-g'(z)=—1 
‘Y=—cota-h’ (z)-iMyg(z)/cosa=—i(k —x, ) 
From the last two equations it is obvious that: 
g(2)=ztana 
h(z)==ztana (x, +ztano/2) 
Assuming a slowly oscillating airfoil, the solution to the PDE has the form: 
9G y) -g() *k[h (2): Myg(z)/cosa ] 
When z = x—-ycoto. 


In order to calculate the pressure distribution on the upper surface of the lower airfoil 


in the zone I the following relation should be used: 

cp =2[ıkp +o] 

Substituting the known functions @ and 0; 

(p=ztana +k{-iztana(x,+ztan’o/2)-(Myztana/ cosa)} 
p,=d/dx[ztana: +k[—iztana(x, +ztan 0/2) GMyztana/cosa))]= 
= tanoiktana[x,+ztan“a—(Myz/cosa)]= 

= tanoik[tano.(x,+ztan 0) +(y/sina-coso) (sina/cosal) ]= | 


=tanoik[tano(x,+ztan/o)+ (y/cos/a)]= 
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=tana—ik[tana(x, +( x-ycotQ)tan at (y/ cos o] = 

=tana-ik[tanau(x, + x tan a.—ytanoct (y/ cos o] = 

=tanoik[tano(x,+ x tan^a—ycota] 

Therefore, the pressure distribution is given by the following relationship: 


cp--2Ítana--ik[tana (x, xtan^Q—yco ta] -ik[ztanact k (-iztana (x, ztan o 2)-AMyztana/ 


cosa}}} 


Taking into consideration that higher order terms of k are neglected the pressure 


distribution is given by the following relationship: 

cp=2tana{-1+ik[x,+ x (tan^a—1)*-2ycota]) 

So at y=0 

cp=2tana{-1+ikfx,+ x (tan’a-1)]} 

In the zone II: 

PEYE) =g(z) +k[Hn(2)M(yg(2)/cosa }+g; (21) +k [hy (z;)+iM(ygi (z1)/cosa] 

where: 

2)=x-Atycota 

In order to calculate the functions g(z),h(z) the boundary conditions will be applied. 
The boundary conditions that must be satisfied at y=h are: 


X, =—cota:p’(z)+cota-g,'(z;) = 0 


'P,=—cota-h' (2) +cota-hy (2,)M[g(2)-g1(z1)]/ cosatiMAcotalg (2) Te (2) ]/2cosa= 0 
After algebraic manipulation: 

cp=—2tana {—1+ik[x)+x (tan^o.—1)] +2A—2ycota} 

For the region III: 

Gy) 7g(z)-k[h (2) 3M (yg(2)/ cosa ] ga (22) k[hiz)) * 3M (yg i / cosod * 

222) k[h2 (22) My gz(z» )/ cosa] 

where: 

Z2—X—À-ycota 

The boundary conditions that must be satisfied at y=0 are: 

X, =- cotocg'( z) +cota-gy'( 21) —cota-go'( z2)= —1 
‘P,=—cota-h’(z)+cota-h,'(z;)—cota-h,’ (z2)-iM[g(z)—21(21) + g2(22)]/cosa = — 1(x —xo) 
After several manipulations: 

cp=-2tana {3+ik[-3x,+ (3x-2A) (tan’a-1)]-2A/cos’a} 


Using the results from the three zones, the pressure jump across the lower blade can be 


calculated from the following relationship 
Acp(x,0)=cp(x,0")- cp(x,07) 


In order to calculate the cp(x,0 ) previous results will be used (two airfoil oscillating case) for 


the lower surface of the store airfoil(store): 
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cp72tana (1-ik[x, x (tan^a.—1)]] 


SoAcp(x,0)=cp(x,0")— cp&,0) 


Two integration intervals can be distinguished in the lower airfoil taking into consideration that 
0.5<A<1: 

1) OSx<A 

2) ASx<2A 


In the first integration interval is: 
Acp,(x,0)=cp(,0°)— ep(x,0 )= 

=-2tana{ I-ik[tanaı(x,+ x(tan’a-1)}= 
=—2tana{ 1—ikftano(x,+ x (tan^a.—1)) and y=0 
Acp,(x,0)=4tana.{ 1-ik(x,+ x (tan^a—1)] 


The latest expression is a complex number with real and imaginary part. In order to calculate 
the pitch-damping coefficient only the imaginary part of Acp is needed. Therefore: 


ImAcp,(x,0)=4tana {—x,+ x (tan^a.—1)) 
In the second integration interval: 
ImAcpy(x0) —-2iktana [-4x;* (4x-2A)(1—t2n^0))-2A /cos o] 
The lift and moment coefficient can be calculated from: 
C, zl 'Acpdx 
=o| "Acp(x-x,)dx 


By definition: 


Gl Ce tik Cy] 
Cue7 Cyert Cue 


For the purpose of the analysis, only the pitch-damping coefficient is required. Therefore, in 
order to facilitate the calculations only the necessary imaginary part of the Cy wil be : 
calculated. 


Cue = Cyr Cyn =(1/ik) of “ImAcp,(x,0) (x-x,)dxtal ImAcpy(x,0) (x-x,)dx= 
o —4tana (x, x (tan "o1 (x—x;)dx-AÍ" 2iktanar[-4xot (4x-2A)(1—tan^0)-2A /cos ^o] (x-x;)dx-— 
--2tna (x; (4-24)-x8A^-4A t 2- (A^-2) (tan’a-1))+2A’-2A-2(2-A))/3 (tan’a-1)] 


The above result is in agreement with the previous results derived by M. F. Platzer and 
H.G. Chalkey but the last one has been derived using the elementary theory. 

Therefore the elementary theory that has been developed predicts the above result as a 
special case for B=0 and 5=180" and represents for that case the pitch damping coefficient of 
an oscillating airfoil mounted close to a stationary airfoil or wall at a distance h/2. It is also the 
same result of the analysis of supersonic flow past a slowly two-dimensional airfoil in a wind 


tunnel with porous wall as presented by M.F. Platzer, 1971. 
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The second general case that will be examined in the thesis is shown in fig. 2. In that 


case, the Mach cone is wide enough to include the wing airfoil. 


Shock 
y Soc SUBSONIC 
LOCUS 

Wing Airfoil 

Reflected 

Shock Waves 

B+3A 
Fig 


Using Sauer’s solution as previously and dividing the flowfield in zones as depicted in 


Fig.7 the perturbation potential @ in zone I is given: 
Pyk) -g(z) tk[h(z)* 1Myg(z)/coso] 

Where: 

z=x-A-B tycota aty=h 

The boundary conditions that must be satisfied at y=h are: 


X, =cota-g’(z) = ES 
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Y, = cota-h’( z)+iMg(z)/cosatiMy g'(z)/ cosa. =~ie? (x-x9-B) 
After integration and algebraic manipulations: 


g(z) =—ztana e 


h(z) =iztanc-e fx, +ztan?0/2+Mh/ COSO) 


Where again the condition has been imposed that the functions g(z) and h(z) are zero for 


negative and zero arguments. 


Ô 


$,——tanoce +iktanor-e | (xp+z,tan°a+Mh/cosa)—My/cosa} 


Ô 


= -tanq e A B o a aNG cosa)-My/cosa} 


=tanoce iktan are { (xo+(x—B)tan a—A tan’atycotath/sina-cosa) 


y/sina-cosa}= 


=-tancre-+iktanar-e? {(x+(x-B)tanart A-ycota} 
Because h=A/cota , tana—1/cosa-sina=cota 


This gives the following result for the pressure distributions on the lower surface of the upper 
airfoil (wing airfoil) in zone I: 


cp =-2{ıkp +9, ] 


cp=2tana-e—2iktana-e® { (xo (x—B) (tan°a—-1)+2A—2ycota} = 


cp= 2tana: e {1—ik { (xp+(x—B) (tan a-1)+2A-2ycota} 
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cp= Zaia A { (xo (xB) (tan2a-1)} 


In order to calculate the pressure on the upper surface of the wing airfoil the same approach 
will be followed and the result is: 


cp- 2tana-e? (—1+ik((xo+(x—B)(tan2o—1)) 


To calculate the Acp of the upper airfoil (wing airfoil) in the interval of zone I we have to 
calculate: 


Acp(xh)=cp(<h )—cp(x;h `) 

Therefore: 

Acp= 4tana-e® {1-1k {(&o+(x-B) (tan^a.-1)) 

Following the same arguments as previously for the zone II: 


Z|7-x-A-ycota. -B 


81 (21) =z,tana-(1-e) 
The second boundary condition gives: 


—cota-h,'(z)+cota-h’(z)iM[-g(z)+g; (z;)]/cosa=-1(x -xo-B-A) 


—cota-h,' (z)*icota-tanae («,+ztan?o.-4+Mh /cos)-iMz,tana/cosat i(x —xy-B—A)=0 
Using the following trigonometric identity: 


ztan’a—z/cos*a=—-z we have: 
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h,(zı) =iz,tana { (xpt z,tan 0/2) Cie e)+Mh ° so) . 


In order to calculate the pressure distnbution on the upper surface of the lower airfoil in the 


zone II the following relation should be used: 
cp =-2[ikQ +0, ] 


Substituting the known functions ( and 9, 


=z,tana- (1-2?) —ztancr-e +k fiz tance {xo+z,tan’o/2) C +e”)+Mhe”/ cosa} rn 


ztan’a/2 +Mh/coso)} 


o,=d/dxfzıtana-(1 pan zytan 0/2) (-1 poco cosa}+ iztana 


go (xou ztanQ/2 +Mh/ cosa).}]= 


=tanor (1-6 0)-tan0c-e +k fitancr{ (xo tz,tan20) (1 +e)+Mhe°/cosa} + en ztan O, 


+Mh/ cosa)} 


= tana-(1 ~2e)+kitanal { (xot z,tan a) ei +e) +Mhe®/ cosa} +e" Got ztana+Mh/ cosQ)] 


Taking into consideration that y=0 we have: 
= i : 2 i iÓ 
=tancr{ (t—2e')+kil (pt 2,tan7a) (-1+2¢)+2Mhe”/cosal} 


cp =—2[ikp +, ]=—-2[tana: (1 ~2e) + kitancrf (xot zıtan d- 2,1)(-1 +2e%)+ 2Mhe®/ cosa}] 


The last relation gives the cp on the upper surface of the lower airfoil (store) in the zone II. In 
the zone III a similar approach will be used to calculate the h(z2) and g2(z2). The variable z2, 
taking into consideration the fig. 2, is given by the following relationship: 

z27x-2A tycota -B | 
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The first boundary condition that must be satisfied at y=h is: 
(oz) =23tan0r-(1-e) 

The second boundary condition gives: 
-cota-h,'(z,)+cota-h'(z)+cota-hz'(z)+1M[g(z)-g:(z1ı)+g2(z2))/cosat4Mhcotafg'’(z)+g,'(z,) 


+9, (2,)]/cosa=-i(&-xo-B-2A) E) 


cota-h; (z,) cola anal) Nao cosaHMh (230% en Bee 
=0 


cota-h; (Z,) = j(x +2 tan 0) +e02,/ cos’a— iMh(2—3e)/ COSQ-1(Z2 —Xo) ES 
h2(z;)7—1iz,tana mana 2e) iMh tana z, (2-36%)/ cosa 
In a similar manner, as previously, the cp can be calculated: 


=z, tana:(1 oe t (1 ES fiz,tana (xy +z,tan“0/2)(1 +e”) +Mhe”/cosa} 
agree (x t ztan“o/2+Mh/ cosa) titana(x,+z,tan’a/2)(1 ¿O ¡Mhtanoz (236%) /cosa 
-iMy[g(2-8161)*:2]/ cosa] 

» ó $ 2 Ó) iÓ iÓ 2 
,=tano-(2-3e )+kitana[(xo+zitan &)(—-1+e )+Mhe /cosa}+e (x,tztana+Mh/cosa)+ 


(xy zan o + 1)-MA (246% /cosa] 


(),=tana: (2-36) +kitanol (xt ztan a) e_wth(2-6e®) /cosa] 


In order to calculate the pressure distnbution on the upper surface of the lower airfoil in the 
zone II the following relation should be used: 


cp=-2[ikp +0, ] 


Substituting the known functions @ and @, 
cp=—2iktana.-[(2-3e)+2, 2-30) + (xt zstan a) -Mh(2-6e)/coso] 


=-2iktana-[(2-3 e +2z2(1 To t(xot (z,-1)tan?o)e®-Mh 2-6.) / cosa] 
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V. METHOD OF CHARACTERISTICS | 


In the method of characteristics, it is desirable to define the coordinate system in such 
a way that all-possible discontinuities could occur across. Along this coordinate system, the 
equations of motion of the flow field can then be treated as ordinary differential equations that 
are solvable by classical or numerical techniques. 

To obtain the equations in this coordinate system (the characteristic directions) in the 


(xj) plane, the equations of motion are wratten in terms of the arbitrary intersecting 


coordinates: 

ç = ¿(x,y) (V-1) 
And 

7 - "x, y) (V-2) 


If the first derivatives of the dependent variables, uy, p with respect to & are made 
indeterminate across lines of T] — constant, and the first derivatives of the dependent variables 
with respect to T] are made indeterminate across lines of Š= constant, then any possible 
discontinuities in the first derivatives will occur across these lines. These lines are then the 
characteristics and their equations are obtained in the evaluation of the indeterminacies. 
Consider first, two-dimensional, steady flow. The governing equations of motion are the 


continuity equation: 


alpu) (pv) _ 
lE Ó (V-3) 


The Euler equation becomes: 
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u——+y—+——=0 (V-4) 


Ox Oy po 
er, oo, (V-5) 


ee (V-6) 
Ox Oy 

Along a streamline: 

aE yy 2 fu 2 +2 (V-7 


Ou av ll-@ à 
de sj Sp ens - 


In terms of the new coordinates: 


= Fr -— 1 
poz us + PÉ Ve + Le, +v b, = -pT),u, — pTyv, -— m. +vn, bo, 
C C 


(V-9) 
In similar fashion the Euler equations become: 
Phi. vé, hu 6p; 7 op, *vn, b, - rn. (V-10) 
pe, «ve, o. Ep, - - pla, *vn, b, - p, (V-11) 
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The last three equations form a system of three equations in uz, Vs, Ps- Solving for pz 


by Cramers rule and taking into consideration that pz is indeterminate across T] both 


determinants in numerator and denominator should vanish. Therefore, from the denominator 


determinant we have: 


[£ +y ]- [Ga —c )6 2 *2vu££, «(Y - e )62]=0 (V-12) 


The solution to the above equation gives the equations of all three characteristics in the 


physical (x,y) plane: 

dy v 

— =—= fan -13 
ang (V-13) 
d 4 

E = tan(¢ +a) (V-14) 


Where € is the angle the streamlines makes with x-axis and æ is the Mach angle. The first 
equation describes a streamline while the last one describes left and right running Mach lines. 


The compatibility equation for u,v, palongé and7 characteristics can be obtained by setting 


the numerator of pz equal to zero: 
O 
K rs ESSA = 0 (V-15) 


The above relation must be satisfied along characteristics (€). 
According to Teipel (1962) the unsteady flow over a flat plate can be treated in a similar 
fashion. In order to facilitate the calculations the continuity and Euler equations can be 


rewritten using the local sonic velocity: 


cs 


o 

l 

= 
O o 


(V-19 


59 


So canceling higher order terms the continuity equation becomes: 


ll ul, CC e ss, (V-17) 
y-lot y—-1 @ Ox Ox 


In a similar, fashion the Euler equations become: 
ðu 2 ðc ðu 


— +——c, —+u,— =0 (V-18) 
Ot y-1 x Ox 


e ie ne ee, (V-19) 
Ot y-l àx Oy 


A system of three differential equations has been formed for three unknowns u, v.c 


To simplify the calculations with the assumption of harmonic oscillation we introduce the 


amplitude functions: 





Hoy > (V-20) 
UL 
V(x, ye = I (V-21) 
E u 
C(x, y)e™ = | omc, 





— vo (V-22) 


Therefore, the continuity and Euler equations become: 


IVE E MEC rc zo (V-25) 
Es ox 
2 ikU=0 ` (V-24) 


L, | EC ikV = 0 (V-25) 


"AMA D 


With the above variables the compatibility relations for the above system of equations are: 








Es. I (a?c-U)- 
E) dr ;— rc v) 0 (V-26) 
ee A E : 
E (E) el ;— rc v) (V-27) 
BUN (ac) | | 
O E) +ikU =0 (V 28) 


Taking into consideration the irrotationality: 


= 0 (V-29) 


The non-dimensional system of equations can be wntten as: 


I. S + M — ue +ikM’C = 0 (V-30) 
E Ox 
OU OC 
— 4 — 4 ikU =0 -31 
Ox Ox (39 
E (V-32) 
oy Ox 


Transforming the last system of equations to the new coordinate system (mn): 


QU OV , M’ 
(82) (8^) La E uno vas 


Ox J; Ox J; 
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(2) (GL i M (v -c)-0 (V-34) 
n n 


Ox Ox ME “El 
(Z) (E +ikU =0 (V-35) 
Ox sir Ox str 


The last system of equations can be solved using finite differences. The values of U,V 
and C on the upper and lower surfaces of the airfoils can be calculated using the results of the 
Elementary Theory (Chapter IIT). The following relations apply: 


U 7 gx, V - tanagy,C = -[px + ikọ] (V-36) 


Using the results of chapter III the following results apply to zone I: 





Ur = tana] - ik(x, 4 xtan? a + ycota)| (V-37) 
M 
Vr una - 1+ik(x — x+ y(cota + »| (V-38) 
y cosa 


C; = tanal- 1+ ik(xo 4 x (imm q — 1) + 2ycot oy (V-39) 


The last three equations represent the values of the parameters on the upper surface of the 
lower airfoil in the interval of the zone I. The last equation gives the pressure distribution and 
it is half the value of the actual pressure distribution coefficient. Following a similar approach 


the following results apply in the zone II* for the pressure coefficient: 


C =2{2tanae?? — tana +iktana[x, + x(tan? a —1)]- (V-40) 
2ike? tana[x, ^ (tan? a — 1(x - B) + 2d cota]} | 
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VI. FLUTTER ANALYSIS 


The lift and moment are given from the following relations: 


1 
pe ; pU?’ |C d (VI-1) 
0 
Iw oe 
M =~ pU [C,G&-x)dx (VI-2) 
0 
And 
2 (00 Op | 
Cp=-—| + — 3 
£ ye | Ot 3 
While 
1 
c, = | Ac dx (VI-4) 
0 
1 
Cy = Í (x — x,)Ac,dx (VI-5) 
0 
And defining: 
c, = 0, [c,,, + ikc >] - (VI-6) 


The last relations have been used in the elementary theory as shown in chapter III, and 
are a result of the following analysis. To compute the non dimensional lift and moment acting 
on the store airfoil the non-dimensional pressure coefficient P as defined in problem 
formulation last equation will be utilized and consequently the lift and moment can be written 


as follows: 
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L= Í {P(x,0*) - P(x,0° ) ix (VI-7) 


ME | {P(x,0° )- P(x,07 a — X, )dx (VI-8) 
0 

And: 

P(x, y) = y.M*C(x, y) (MIES) 

Where: 

Kae (VI-10) 


e 


Although in the thesis only pitching, oscillation is considered the more general case 
include plunging oscillation also. In that case, Garrick and Rubinow’s (1946) method of 


expressing lift and pitching moment is used: 


ES Sp Pe (L, +iL,)+0,(L, +iL,)] SEND 
C 


29 ; | 
M= EL ukit o cM, + iM.) + 69(M3 +iM4)] (VI-12) 
c 


Where: 


=> 


qe aade = 2b 
u 


e 


For a single degree of freedom as in that case the flutter analysis is much easier, and the 


equation of motion is: 


1,0 +(1+ig)C,0 = M, (VI-13) 


The sum of the moments of inertia about the elastic axis 1s: 

M, =-1,0 (VI-14) 
In addition, the elastic restoring moment is: 

M, =-(+ig)C,0 (VI-15) 


Moreover, the aerodynamic pitching moment about the elastic axis 1s: 


M, = -2 p. cu 2 k2e'20 (M, + iM,,) (VI-16) 
Since: 

0 =0,e“ (VI-18) 
Ó=-a?0,e” (VI-19) 


So the equation of motion becomes: 

| 2 
e'” [-1,0”@, +(1+ig)C,0, = -4p c u, k’e'*0 (M, +iM,) (VI-20) 
Separating the equation in real and imaginary parts gives: 


- Io! * C, + pit M, = 0 (VI-21) 


we 
gC, + a p.c a M, =0 (VI-22) 


Using the results of Garrick and Rubinow (1946): 


Q,X-ur, +M, =0 (VI-23) 
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M,+22,X =0 (VI-24) ^ 


In order to solve the single degree of freedom problem it is necessary to obtain values 
of M, as a function of k untl the second equation is satisfied. The value of M, for that k is 


then used 1n the first equation to determine if X is reasonable: 





X=1-2>0 (VI-25) 


The non-dimensional flutter frequency is: 


DER 





=— (VI-2 
Do a o 
Accordingly the speed is: 
er l 
— = — 1:2 
6o, kd X Mes 


Aerodynamic instability (flutter) occurs when M,S O while the L, M coefficients are given by 


the following: 
2 ] } 
L=L, +iL, =-G1 [C(x,0*)ax - [C(x,0" ax} (VI-28) 
0 0 
L 


M - M,*iM, -- [C(x,0° ns S D 


And defining the following quantity, the non dimensional pitching moment, the analysis can 
be greatly facilitate: 


OO Ice rire, les (VI-30) 
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(VE31) 


(VI-32) 





f 
D 





VII. CONCLUSIONS 


In this thesis, Platzer’s (1973) approach to analyze the aeroelastic stability of flat-plate 
cascade and of airfoils in wind tunnels was extended to the case of two airfoils in close 
proximity to each other. The airfoils could oscillate in pitch with an arbitrary phase angle 
relative to each other. The analysis is valid for supersonic flow and is based on the small 
perturbation theory. 

Two.cases needed to be distinguished, namely airfoils with supersonic leading edge 
locus versus subsonic leading edge locus. In the former case, the leading edge of the upper 
airfoils is located upstream of the Mach wedge generated by the lower airfoil. In the latter case, 
the upper airfoil is within the Mach wedge of the lower airfoil 

For the case of supersonic leading edge locus a general formula could be derived for the 
pitch damping of the lower (store) airfoil which gives the dependence on Mach number, pitch 
axis location, phase angle, and proximity of the two airfoils. For the special case of out-of- 
phase oscillation between the two airfoils the result agrees with that derived by Platzer and 
Chalkey (1972), equations (37). The results showed that interference has a strongly 
destabilizing effect. In addition, for the same special case Platzer’s (1973) equation (4.2) could 
be verified. 

Furthermore, the analysis of the subsonic leading edge problem was begun and the 
pressure distribution on both sides of the upper foil was derived. In addition, a method of 


characteristics procedure was described which could be used to compare results. 
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